Aging Effects of an Elastic String Diffusing in a Disordered Media 
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The aging effects of a 'diffusing glass', a single elastic string diffusing in a two-dimensional disordered 
media assisted by thermal noise are studied by Monte Carlo (MC) simulations. We find for the 
first time convincing numerical evidence of non-trivial aging effects both in the linear response 
against transverse force and the associated correlation function. Our results retain some important 
predictions of the conventional dynamical mean-field theory (MFT) but also reveal remarkable 
differences from it. 

PACS numbers: 74.60. Ge, 02.50.Ey, 75.50.Lk 

Elastic objects diffusing in disordered media assisted 
by thermal noise, which we call simply as 'diffus- 
ing glasses', are expected to remain persistently non- 
stationary due to the complicated competitions between 
the elasticity and disorder. An example will be a single 
vortex line diffusing in a type-II super-conductor with 
randomly distributed pinning centers. A dynamical MFT 
of elastic manifolds in random media pi] has been devel- 
oped which proposes a self-consistent picture for such dy- 
namics and predicts existence of aging effects similar to 
those found in spin-glasses. However, the MFT is exact 
only in infinite dimensional embedding space. Compared 
with many driven dynamics || of the same elastic ob- 
jects, much less is known about such diffusive dynamics 
in realistic dimensional spaces. 

In order to go beyond the MFT, we study in this let- 
ter the dynamics of a simplest 'diffusing glass', namely 
directed polymer in random media (DPRM) Q on a two 
dimensional square lattice by extensive MC simulations. 
The previous MC simulations [j|]6) have found encour- 
aging signatures of non-stationary dynamics under zero 
driving force but many open questions remained. 

We find that the k = component of current (veloc- 
ity) correlation function and linear conductivity exhibit 
novel scaling behaviors which reveal 'one step' violations 
of the time translational invariance (TTI) and the Green- 
Kubo (GK) formula or the fluctuation dissipation theo- 
rem (FDT) for 'currents' ]7j]. On the other hand, finite k 
components equilibrate in the time Tfc oc k~ z ( T \ 

Model- Our model is a lattice string on a two dimen- 
sional square lattice of size L in the longitudinal direc- 
tion z and M in the transverse direction x. The string 
is directed in z direction and its segments are labeled as 
x(z) (z = 1, 2, • • • , L). They obey the RSOS (Restricted 
Solids on Solid) condition, i.e. \x(z) — x(z — 1)| =0, ±1. 
Random potential V(z,x) is defined on each lattice site 
(z,x) which takes a random value from the uniform dis- 
tribution between -1 and 1. This system is known to be 
in the glassy phase at all finite temperatures. 0] 

The dynamics is introduced by a heat-bath type MC 
method |3]|| which ensures that the microscopic motions 
of the kinks and anti-kinks are thermalized with the heat- 



bath at a temperature Tb a th- In one Monte Carlo Step 
(MCS), the whole configuration is sweeped once. The 
energy of the system at time t is given by 

L 

E(t)[V,h] = ^T f V{z,x(z,i)) - h(z,t)x(z,i), (1) 

z=l 

where h(z,t) is the transverse force. As in the case of 
static properties pj , we assume that the dynamics of the 
present lattice model and the continuous model studied 
in the MFT have same asymptotic scaling properties. 

We have performed simulations similar to the zero field 
cooling (ZFC) experiments of spin-glass systems. 0]. At 
first, a certain waiting time t 2 is elapsed under zero field 
starting from an out-of equilibrium initial configuration. 
Then a constant filed h(z') is applied at z' afterwards. 

If linear response holds, the induced current at z, i. e. 
temporal transverse velocity of the segment x(z), to be 
measured at time t\{> t 2 ) can be written as, 



SJz-z'{ti,t 2 ) = / dt'a z 
Jt 2 



'(h,t')h(z'), 



(2) 



where cr z - z >(t,t') (t > t') is a time-dependent conduc- 
tivity. Instead of the currents, we measure the linear 
susceptibility (induced displacement divided by h), 
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dt 



5J Z - Z ,(t,t 2 ) 
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(3) 



Let us here introduce a generalized Green-Kubo (GK) 
formula. 
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< J z (ti)Jz>(t 2 ) > e(h -ta), 

(4) 



where < J z {ti)J z >(t 2 ) >= d tl d t2 < x(z, ti)x(z', t 2 ) > is 
the current correlation function. The usual GK formula 
corresponds to the case with the 'FDT ratio' Y = 1. 
Hereafter, the bracket < • ■ ■ > means the average over 
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samples: different realizations of initial configurations, 
thermal histories (MC runs) and random potentials. 

By integrating the current correlation function over the 
two time variables we obtain, 



B, 



'(*i,ta) 



dt / dt' < J z {t)J z ,{t') > 

t 2 Jt 2 

< [x(z, t\) - x(z, t 2 )] [x(z', ti) - x^z', t 2 )] > 



Combining with (||), the integral violation of the FDT 
can be defined as, 



>(h,t 2 ) 



B z 



>(h,t 2 )/2 - T ha . thX z-z'(ti,t 2 ). (6) 



We define the Fourier transform along z axis as 
Xk(ti,t 2 ) = Y,r=o cos ( kr )Xr(ti,t 2 ) and B k {tx,t 2 ) = 
Sr=o cos(kr)B r (ti, t 2 ) where k = mr/L with (n = 
0, 1, • ■ ■ , L). Note that B k= o is proportional to the mean- 
squared displacement of the center of mass. 

In the following, we begin with the crossover between 
the linear and non-linear response. Then, we discuss ag- 
ing effects in linear responses and correlation functions. 

Crossover between linear and non-linear response.- As 
often emphasized, the non-linear effects due to the collec- 
tive creep || should be dominant asymptotically. How- 
ever one can systematically expel them out of a given 
time window by decreasing h because the characteristic 
time to create the nucleus of the creep diverges rapidly as 
h is decreased by the well known formula lnr crccp oc hr^ 
with the glassy exponent /i = 1/4 in two dimension. We 
demonstrate it in the following. 

We have performed simulations of system size L = 500 
M = 520, at temperature Tbath = 0.4 under uniform 
transverse fields. The periodic boundary condition is im- 
posed on both x and z directions and straight lines are 
chosen as initial configurations. 

In Fig. 1, we show data of the susceptibility X k=o(j, 0) 
at different strength of fields h in a double logarith- 
mic plot. There are crossovers from the linear response 
curve where data of different h merge with each other to 
the non-linear branches which depart from the common 
curve. The non-linear branch grow linearly with r by 
different velocities v(h). In the inset of Fig. 1 we show 
the effective velocity v(h) determined by fitting the non- 
linear branches to the form 1ix(t, 0) = v(h)r + c(h) where 
c(h) is a parameter. The result is consistent with the ex- 
pected behavior hiv(h) oc — lnr crc0 p. More details of the 
crossover phenomena will be reported elsewhere. 

Aging Effects - Let us begin discussions on aging effects 
in the linear susceptivility \ and integrated correlation 
function B at zero wavenumber k = 0. 

In Fig. 2 we show the data of Xk=o and Bk=o of differ- 
ent waiting times t 2 against the time difference r = t\—t 2 
in double logarithmic plots. The system size is L = 500 
and temperature is Tbath = 0.4. We have checked that 
there are no finite size effects up to 10 6 MCS within the 



statistical accuracy by simulating larger systems. We 
have chosen weak enough uniform field h = 0.005 in or- 
der to measure linear responses up to 10 6 MCS. 
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FIG. 1. Crossover between linear and non-linear responses. 
The strength of the uniform field is varied as h — 0.64 ,0.32, 
0.16, 0.08, 0.04, 0.02, 0.01, 0.005, 0.0025 from the top curve to 
the bottom. The error bars are of the size of the symbols. The 
average is took over 100 samples. (1000 samples for h — 0.005 
and 0.0025) The inset is the velocity v(h) versus /i -1//4 at 
h = 0.64 ,0.32, 0.16, 0.08, 0.04. 

Except the rapid growth at r < 10 2 MCS, which is 
discarded in the following as a short time transient be- 
havior, the general feature is the following: each curve of 
a given waiting time t 2 follows initially the lower 'quasi- 
cquilibrium' branch then switches over to the upper 'ag- 
ing' branch at around logr ~ log £2- Here the existence 
of aging effects is evident. Both \ and B violates TTI: 
they are not functions of the time difference r = t\ — t 2 
alone but explicitly depends also on the waiting time t 2 . 
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FIG. 2. rbttthXfc=a(*i)*2) and B k =o(ti,U)/2 of waiting 
times t a = 10, 10 2 , 10 3 , 10 4 , 10 s , 10 6 from the left curve 
to the right. Only lines connecting the data points at 
r = ti — 12 ~ nx 10 p (where n = 1, 2, • • • , 9 and p = 0, 1, • • • , 5) 
are shown for graphical convenience. The average is took over 
1050 samples. Typical sizes of the error bars are those of the 
symbols on the curve of t% = 10 6 (bottom). 
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FIG. 3. Scaling plot of the susceptibility and the 
mean-squared displacement, (lower and upper curves respec- 
tively) The inset is the scaling parameter xea(^2) vs ti. We 
use the data of r > 10 2 (MCS) at t 2 = 10 3 ,10 4 ,10 5 ,10 6 (MCS). 
All the broken lines are algebraic law fits of the same expo- 
nent 2/z(Tb a th = 0.4) = 0.42 with different amplitudes, from 
which we obtain cijc\ ~ 4.0 and y ~ 0.5. 
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FIG. 4. Parametric plot of the scaling functions x(i) and 
B(t) shown in Fig. 3. The solid straight line represents the 
integrated FDT. The dashed line is a fit by a straight line 
with slope y = 0.54 to the part where the FDT is violated. 

The one step structures can be well described by the 
following simple scaling ansatz: 



IbathXfc=o(tl)*2) 

B k=0 (t u t 2 )/2 



S x (ti/t 2 )(*i 
S B (ti/t 2 )(ti 



.t 2 )W). 



(7) 



Here 5 x (t), 5s(£) are step-like functions in the scaled 
time variable t = ti/tz. They take a) the same value 
ci in the 'quasi-equilibrium time domain' i < 1 and b) 
different values yc2{> c±) and c 2 (> ci) respectively in 
the 'aging time domain' ( > 1. The ratio y is smaller 
than 1. The latter means that the integral violation of 



FDT (|6|) exhists : / > in the 'aging time domain'. 

In the absence of the random potential, the present 
model is equivalent to the free Gaussian field model @] 
in which z — 2 : the diffusion is 'normal'. Further more 
S x and 5b are constants (mobility or diffusion constant) 
and equal. Thus TTI and the FDT holds completely. 

In order to test the scaling ansatz (]?]), we consider 
equivalent formulas: T ba thXfe=o(*i. h) - XEA(h)x(t) and 
B*=o(*i,*a)/2 ^ XBA(t2)B(i) where xea(^) = t 2 2 /z(T) 
scales the 'height of step'. The scaling functions read as 
X(t) = S x {t)\t-l\ 2 ' z ^ and B(t) = 5s(?)|?-l| 2/z(T) . We 
assume that i — t\jt 2 is the correct scaling variable and 
treat xea(^) as a scaling parameter || to be determined 
for each t 2 . The resultant master curves and the scaling 
parameters xea(*2) are shown in Fig. 3. In Fig. 4, we 
show a parametric plot of the master curves which clearly 
shows 'one step' violation of the FDT. 

From (Q), we can deduce scaling forms for the linear 
conductivity and correlation function of the current at 
k = 0. Except around the step \og(ti/t 2 — 1) ~ 0, we find 
power laws multiplied by the 'one step' functions, 

< J(h)J(t 2 ) > * s B {ti/t 3 ){h - h)-^ 1 - 1 '^. (8) 

In the last equations we have used o~{ti,t 2 ) = 
-d tl d t 2X{h,h), < J{h)J{t 2 ) >= -d tl d t2 B(t u t 2 )/2 
and omitted the common pre-factor — 2/z(T)(l— 2/z(T)). 
So both noise and response violate TTI and the FDT. 
The 'FDT ratio' Y(ti,h) (@) shows 'one step' variation 
: Y = 1 in the the 'quasi-equilibrium time domain' and 
Y = y in the 'aging time domain'. 

Let us now compare our results with the conventional 
picture based on the MFT H . There are now increasing 
number of numerical studies on different glassy systems 
such as a spin-glass model JTo) , Lennard- Jones Glasses 
jn]], and coarsening in spin-systems fi^l , which appear 
to roughly support the picture based on similar dynamical 
MFTs 131. However we find that the crucial concept of 
the MFTs Jl],[l3| called as 'correlation scales' cannot be 
applied to our system as we discuss below. 

The MFT assumes that TTI and the FDT hold in the 
quasi-equilibrium scale B < B^a but not in the aging 
scale B > Be a- The border line between the two scales 
Sea is a well-defined constant. Consequently, the cor- 
relation and response functions becomes the sum of the 
contributions from the quasi-equilibrium and aging scale. 

However, the multiplicative scaling form (0) certainly 
disagree with the MFT. The 'height of step' Xea(*2), 
which corresponds to the Bea of the MFT, is not at all a 
constant but apparently increases with time t 2 as shown 
in the inset of Fig. 3 (similar effect has been noticed in 
the Sinai model H])- Thus our results reveal consid- 
erable differences from the MFT. However, surprisingly, 
the 'one step' variation of the FDT ratio predicted by the 
MFT is recovered in an unexpected way. 



3 



Propagation of response - Finally we discuss the prop- 
agation of responses along z axis. To this end, we have 
performed another set of simulations with point fields 
h(z) = h po intS(z — L/2) which pull the center of the string 
z = L/2. We have obtained Xr(h>h) by measuring the 
induced displacement of segments at various distance r 
from the center. 

In Fig. 5. we show a set of data in a scaling plot 
in order to test the scaling law which holds in the pure 
system @, 

Xr(h,t 2 ) = Xr=0(tl,t 2 )H{r/ X r=0(t 1 ,t 2 )). (9) 

Here Xr=o(ti,t 2 ) plays the role very similar to the 'do- 
main size' in the coarsening systems fl5| . In the scaling 
analysis, we have treated each Xr=o(ti,t 2 ) as a scaling 
parameter to be determined for each (ti,t 2 )- 
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FIG. 5. Scaling plot of the profile Xr{U,t 2 ) at t 2 = 10 3 , 
10 4 , 10 5 , 10 6 and r = 10 3 , 10 4 , 10 5 , 10 6 . The system size is 
L = 500 and temperature is Tbath = 0.4 The average is took 
over 11520 samples. The non-linear effect is much weaker than 
in the case of uniform fields and we have used safely point field 
of strength h po i n t = 0.1 to investigate linear-response up to 
10 6 MCS. 

The resultant master curve H(f) drops linearly H(f) = 
1 — sf at f < 2 and has a Gaussian tail at f > 10. Surpris- 
ingly, it is virtually indistinguishable from that obtained 
by simulating a corresponding pure system within our 
numerical accuracies. 

In the Fourier space, the scaling form (g) takes 
the same form as in the MFT |Q, k 2 Xk(ti,t 2 ) = 
F(k 2 Xk=o{h,h)) where Xk=o(*i,*2) = hXr=o(h,t2). 
(Our data satisfy latter within our numerical accura- 
cies.) In the last equations we used h = J °° dyH(y) and 
F(x) = (x/h) J °° dy cos(j/ \J x / h) H (y) . Combining with 
( f7|) we find th at 'domain size' grows as Xr=o{ti,ti)) — 
y/S x (t 1 /t 2 )/h(t 1 - t 2 ) 1/z{T) . In other words, a finite k 
component equilibrates in the time oc fc~^ T ) to the 
equilibrium value x^ K k~ 2 . The latter k~ 2 scaling is 
due to the so called statistical tilt symmetry |l6| . 

The exponent z(Tb a th = 0.4) ~ 4.8 is much larger than 
in the pure model z = 2. We have found that z(T) in- 
creases with decreasing temperature. Though the previ- 



ous work H proposed a logarithmic law for the suscep- 
tibility the present data of increased statistical accuracy 
over enlarged time range fit better to the algebraic law 
m (0). Thus the naive scaling argument based only on 
typical energy barrier || which suggest the logarithmic 
law should be inaccurate. 

To summarize, we performed extensive MC simulations 
and scaling analysis on the aging effects of a 'diffusing 
glass', 2 dimensional DPRM. We found that the corre- 
lation and linear response functions exhibit interesting 
'one-step' scaling phenomena. It will be very interesting 
if such features can be observed experimentally in some 
'diffusing glasses', such as vortex lines in dirty type-II 
super-conductors. 
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